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Dipolar ordering of random two-dimensional spin ensemble 
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We theoretically study the randomly positioned two-dimensional system of interacting magnetic dipoles rep- 
resenting the monolayer arrays of single-domain particles. It is showed the onset of the dipolar in-plane 
superferromagnetic ordering of Ising spins at the surface concentrations of nanoparticles above 0.24. The 
critical concentration of particles with random easy axis arrangement at zero temperature is 0.65. It is 
demonstrated that the ensemble with Ising arrangement of spins is ordered at high temperatures close to a 
particle Curie point. 
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Dipolar ferromagnetism, due to interactions between 
the magnetic moments of rare-earth atoms in crystal lat- 
tices, was first observed at very low temperatures 1 . Re- 
cently, ordering associated with dipole-dipole coupling 
between monodomain particles randomly placed on a 
surface was discovered at high temperatures 2-4 . The 
nanoparticlc in such systems in view of its small size acts 
as a single spin. 

Zhang and Widom 5 applying generalized mean-field 
theory revealed that the three-dimensional random sys- 
tem of spins within infinitely prolated sample can have 
dipolar ordering. However, the dipolar ferromagnetism in 
a moderately elongated spheroid is impossible 6 . From ge- 
ometrical considerations, one can expect in-plane dipolar 
ferromagnetism in the two-dimensional systems. Using a 
point source model which is only valid at the vanishing 
spin concentrations, Meilikhov and Farzetdinova 7 showed 
the absence of the dipolar ordering in the random two- 
dimensional system of Ising dipoles. Nevertheless, it was 
demonstrated in Rcf. 7 employing Monte Carlo modeling 
a pcrcolativc transition to the dipolar ferromagnetism in 
square lattice. Scheinfein et al. 8 studied by the mean- 
field theory and Monte Carlo simulations the ordering 
of random Fe islands on the surface and concluded that 
dipole-dipole interactions are insufficient to explain room 
temperature magnetism observed. 

In this work, by means of generalized mean-field the- 
ory, we show the possibility of in-plane dipolar ordering 
in the two-dimensional sufficiently closely packed system 
of the spins. 

Let us consider the random two-dimensional monolayer 
system of identical single-domain particles. In the first 
approximation, we can suppose the nanoparticles as mag- 
netic dipoles which interact with each other spin. At sur- 
face concentrations of cylindrical nanoparticles c = Ns / S 
(s is the area of the particle, S is the area of the sample, 
N is the number of the particles) as high as 0.6 projection 
H of the net random field of all the dipoles onto selected 



direction satisfies the Gaussian distribution 9 : 
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where r is the two-dimensional radius- vector, r (m) is the 
distribution function over single-domain particle mag- 
netic moments m, whose absolute values are supposed 
to be equal to M s v for all the particles, M s is the satu- 
ration magnetization, v = sh is the volume and h is the 
particle height, s = irr^, tq is the particle radius, 
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Integration over radius r in Eq. 2 begins from 2ro, which 
is the minimal possible distance between two particle cen- 
ters. 

We first consider the in-plane Ising model. The mag- 
netic moments of the particles may have two opposite di- 
rections lying in the plane of the monolayer, Na dipoles 
are oriented along a selected direction and N (3 have op- 
posite alignment, a + (3 = 1, m = a — [3 is an order 
parameter. For Ising model, the distribution over mag- 
netic moments is given by 



r(7) = a5{"f) + 06^ — tt), 



(4) 



where 7 is the angle between m and the selected direc- 
tion. 

The parameters of the Gaussian distribution are read- 
ily obtained: 



#011 



M s vircm hocm 
2rns s 



2„,2 



'II 



16s(2r ) 4 



M2- 



s 

(5) 

The order parameter can be found from the next equa- 
tion, 
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where T is the temperature and ks is the Boltzmann 
constant. 

The most interesting critical parameter for such a sys- 
tem is critical concentration above which the ensemble 
has superfcrromagnetic ordering. Assuming that m is 
small in Eq. 6, we obtain the equation for critical con- 
centration c c : 
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HdH = 1. 



(7) 

This equation can be also used for the calculation of the 
critical temperature at fixed concentration. At zero tem- 
perature, tanh(M s vH /ksT) in Eqs. 6 and 7 may be re- 
placed by sign function sgnH. In this case, it is possible 
to derive the analytic relation for the critical concentra- 
tion: 
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Numerical value for cq does not depend on material and 
is equal to ll7r/128 ps 0.27. As we can see concentration 
for exploiting the Gaussian distribution c « 0.6 exceeds 
this value. In order to refine calculations, we can repre- 
sent the negative logarithm of a characteristic function 
as a series expansion and obtain the distribution over the 
dipolar field after inverse Fourier transform 9 , 
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where moments p n of function r(m) are defined as fol- 
lows: 



/!„ = J ^™r(m)dmdr. 



(10) 



In our computations, we limited the expansion by the 
tenth degree of p. It should be emphasized that odd- 
numbered moments p„ depend linearly on m. Substitut- 
ing Eq. (9) into Eq. (6) and assuming that m is small 
near the critical point, we deduce following equation to 
calculate c c : 
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FIG. 1. Phase diagrams for in-plane Ising dipoles. The areas 
under the curves correspond to the ordered state, T c is the 
Curie temperature of nanoparticle material, Td is the critical 
temperature of dipolar ordering. 



For the in-plane Ising ensemble of spins the moments are: 
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Eq. (11) can be solved numerically. Using the 
parameters of the magnetic nanoparticles observed 
experimentally 2,4 , we calculated that for Co particles 
co ~ 0.2376 and for magnetite Co ~ 0.2374. These values 
differ little from cq obtained using the Gaussian distribu- 
tion. 

The phase diagrams showing conditions for the dipo- 
lar alignment of in-plane Ising spins are illustrated in 
Fig. 1. The curves were calculated by virtue of solving 
Eq. (11). For the numerical simulation of cobalt parti- 
cles, we utilized the dependence of M s on temperature 
proposed in Ref. 10 with Curie temperature T c = 1385 K. 
The magnetite nanoparticles was modeled using expres- 
sion for M S (T) based on the Weiss theory: 
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It is evident from Fig. 1 that temperature of the Ising 
ensemble dipolar ordering is mostly constrained by the 
Curie point of the particles. 
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Further, we consider the two-dimensional ensemble of 
the single-domain particles with random orientation of 
their easy axes. We suppose that the selected direction 
(z-axis) again lies in the plane and the particle magnetic 
moments may acquire two opposite directions. Now, Na 
is the number of spins pointing into half-space with pos- 
itive z and N P dipoles point into another half-space, 
m = a — f3. In this case, the distribution over magnetic 
moments is 



r(j,ip) 



a sin(7) — /3 cos(7) 
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where tp is the azimuthal angle of m, 7 varies within 
range [0, tt/2] and the equation for determining magneti- 
zation becomes 



tt/2 



sin 6 d6 



, M s vH cos# Trr ,„ 
tanh — i- W{H) dH 



(13) 

After the integration over 9 we arrive at following equa- 
tion: 
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The equation for the critical concentration calculation 
with Gaussian as W(H) turns to 
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The parameters of the Gaussian distribution of the sys- 
tem of the particles with the random easy axis orientation 
are 
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At zero temperature we come again to formula (8), 
which gives now cq = 57r/24 w 0.65. This value is ap- 
propriate for making use of the Gaussian distribution of 
the dipolar field. Fig. 2 shows the phase diagrams calcu- 
lated with Eq. 15. We can see from this figure that the 
critical temperatures are below than in the case of paral- 
lel Ising dipoles. This can be attributable to the chaotic 
arrangement of spins and, hence, their fields preventing 
ordering. 

Fig. 3 depicts the order parameters of magnetite parti- 
cle system as a function of temperature calculated using 
Eqs. (6) and (14). Obviously, the degree of magnetic or- 
dering and the critical temperature of random axis spin 
ensemble lie much below than ones of the Ising spin sys- 
tem. The experimentally observed distribution of mag- 
netite nanoparticle magnetic moments was between the 
random axis and Ising arrangements 4 . Hence, our calcu- 
lations do not contradict these experimental data. 

Ferromagnctically ordered two-dimensional nanostruc- 
tures examined in Rcfs. 3,8 have surface concentrations 
c w 0.5 which are below cq of the random axis model. 
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FIG. 2. Phase diagrams for dipoles with random orientation 
of axes. The areas under the curves correspond to the ordered 
state, T c is the Curie temperature of nanoparticle material, 
Td is the critical temperature of dipolar ordering. 
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FIG. 3. Dependence of the order parameter of magnetite 
nanoparticle two-dimensional ensemble on temperature for 
c=0.8. 



Scheinfein et al. 8 pointed out that exchange is responsi- 
ble for the order in their experiments. It seems that two- 
dimensional Co nanodot assemblies fabricated by Fang 
et al. 3 has anisotropy due to growth process. Monte 
Carlo algorithm 11 used in Rcf. 3 also accounts for the 
anisotropy. Thus, these data do not quite correspond 
to the random axis model which has the higher critical 
concentration for dipolar ordering onset. 

It worth mentioning that in the case of spins oriented 
perpendicularly to the plane Hq± > (Rcf. 9 ) so that 
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only antiferromagnetic ordering may exist. Obviously, 
this fact arises from the demagnetizing field of the pla- 
nar dipole system which has maximum in the transverse 
direction. It should be also emphasized that the in-plane 
order in the two-dimensional system has local charac- 
ter limited by a certain size 12 . The experimental results 
show splitting into domains 4 . 

In summary, we theoretically demonstrated the possi- 
bility of dipolar fcrromagnctism in the two-dimensional 
random system of spins. We found the critical concen- 
trations of finite dipoles at which this order sets in. It 
follows from the results of the modeling that the system 
with in-plane anisotropy has dipolar order at high tem- 
peratures. 

The author is grateful to Dmitri Feldman for useful 
discussion. 
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